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TORSIONAL STIFFNESS OF THIN-WALLED SHELLS HAYING REINFORCING CORES AND 
RECTANGULAR, TRIANGULAR, OR DIAMOND CROSS SECTION 1 

By Habvey G. McComb, Jr. 


SUMMARY 

A theoretical investigation has been made oj the Saint- 
Venani torsion of certain composite bars. These bars are 
composed of two materials — one material in the form of a 
thin-walled cylindrical shell and the other material in the form 
of a core which fills the interior of the shell and is bonded to it. 

An approximate boundary-value problem is formulated on 
assumptions similar to those of the theory of torsion of hollow 
thin-walled shells ( Bredt theory). This boundary-value prob- 
lem is solved exactly for a rectangular cross section and approxi- 
mately for slender triangular and slender diamond cross sections. 
Results for the torsional stiffness constants are presented 
graphically. 

INTRODUCTION 

Certain airframe components such as wings, stabilizing 
fins, control surfaces, and helicopter rotor blades have been 
fabricated by employing a high-strength shell bonded to a 
core made of some fight-weight material. The shell is 
formed in the external contour of the component and the 
core fills the interior and acts to stabilize the shell against 
local buckling. Such a structure has been called a “foam- 
filled shell” because the core is often a foamed-plastic 
material. Metal honeycomb and balsa wood have also been 
used for cores. 

A large amount of literature exists on the problem of 
torsion of homogeneous isotropic cylindrical bars, but rela- 
tively little work has been done on the torsion of composite 
cylinders. A few exact solutions to problems in the torsion 
of composite sections are presented in references 1 and 2. 
Solutions for other cross-sectional configurations are desir- 
able, and the methods used in references 1 and 2 do not 
appear to be applicable for sections of the type considered 
in the present report. 

In this report an approximate boundary-value problem 
for the torsion of foam-filled shells, is formulated. The 
fact that the thickness of the outer shell is small relative to 
overall dimensions of the cross section allows an approxima- 
tion similar to that of the Bredt theory for the torsion of 
hollow thin-walled shells with free warping. (See ref. 3, 
pp. 298-302.) For a rectangular cross section, an exact 
solution to the approximate boundary-value problem is 
obtained. For slender triangular and slender diamond cross 


sections, approximate solutions which appear to be reasonably 
accurate are obtained. Finally, the results are compared 
with results based on an elementary concept of the torsional 
stiffness of foam-filled shells. 


SYMBOLS 


A,B 

arbitrary constants 

a,b,c 

cross-sectional dimensions (see fig. 4) 

C 

curve defining boundary of a region 

/(«)*(© 

arbitrary functions 

Qi,Qt 

shear moduli of shell and core materials, 


respectively 

K(x) 

Fourier coefficient (see eq. (21)) 

Ji,J* 

torsional stiffness constants 

fl to 

„ Oib h2 2 

*-~Qit~ G l t 


hinf^2n 

arbitrary constants 

M 

moment on cross section 

m 

index in equations (37) and (61) . 

n 

summation index 

R 

region 

s 

direction tangential to a curve in cross 


section 

t 

thickness of shell wall 

to 

dimension of cross section (see fig. 4) 

U 

total complementary energy per unit 


length 

u=g 2 u 


x,y,z 

coordinates along X-, T-, and Z-axes, 


respectively 

<*«,P 

arbitrary constants 

7 

shear strain 

e—b/a 

e 

angle of twist per unit length 

X, 

eigenvalues 

M = Vl + £S 


V 

direction normal to a curve in cross 


section 

Z,V 

nondimensional coordinates in x- and 


y-directions, respectively 

T 

shear stress 


stress function 

Q 

constant of integration 


i Supersedes NAOA Technical Note 3749 by Harvey G. McComb, Jr., 1956. 
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refer to regions and curves indicated in 
figure 1 

Prime denotes differentiation with respect to the indicated 
variable. 

BASIC EQUATIONS , 

In this section equations of elasticity are established for 
the Saint-V enant torsion of certain cylindrical bars composed 
of two materials. These equations are then particularized 
to the case of foam-filled thin-walled shells. Finally, the 
energy approach to the torsion of foam-filled shells is dis- 
cussed briefly. 

TORSION OF CYLINDRICAL BARS CONSISTING OF TWO MATERIALS 

Consider a long cylindrical bar composed of two. iso tropic 
materials in which one material surrounds the other. A 
cross section of such a bar is shown in figure 1. For torsion 
with free warping, the stresses are given in terms of a stress 
function <f> as follows: 


_2 

Subscripts: 

1,2 




>- 


J 


( 1 ) 


where t=l or 2. 
equation 


Each function <f> { must satisfy Poisson’s 
V 2 4>,=-2 G t 0 (2) 


in its corresponding region R { . 

The boundary conditions can be expressed in terms of the 
stress function by consideration of the components of shear 
stress normal and tangential to a curve in the. cross section. 
These components are, respectively, 



Fia toe 1 . — Cross section of cylindrical bar composed of two materials. 


The condition that the external boundary of the section must 
be free of stress is obtained by integrating the first of equa- 
tions (3) along C\: 


0 i| Ci =Qi= Constant (4) 

The conditions which must be satisfied at the interface C% 
can be seen by referring to figure' 2. The tangential strain 
must be continuous across C\. In terms of the stress func- 
tion, this condition is 


1 < 1 

6i dv c } Qi dv 

Lastly, the shearing-stress component normal to CJ 2 must be 
continuous across C 3 , or 


(5) 

c. 




( 6 ) 


When equation (6) is integrated, the following equation 
results: 

( 7 ) 


The problem is to find stress functions <j > i and which satisfy 
equation (2) in their respective domains and the boundary 
conditions (eqs. (4), (5), and (7)). 

The total moment on the cross section is given by 

f f (T^—Tx.yWRt (8) 

s { 


This equation can be written in terms of ^ as follows: 

<9) 


z 



Figure 2. — Shearing stresses on an element at the interface between 
the materials. Superscripts on symbols correspond to rogions 
indicated in figure 1 . 
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Integrating by parts and making use of equations (4) and (7) 
gives 

M=T2 f (—x dy+y dx)+J2 f f 2 fa dR t (10) 

7=1 t/ Of 4/ »/ 


The stresses and moment on the cross section are independent 
of fl 4 ; therefore, these constants may be chosen arbitrarily. 
For convenience, and Qj are both taken to be zero so that 
the expression for the moment is analogous to that for the 
torsion of a uniform cross section; that is, moment equals 
twice the volume under the 4> di a gram. 


TORSION OF FOAM-FILLED SHELLS 

The equations of elasticity are particularized for the case of 
a cylinder made of a thin-walled shell of one material filled 
with a core of another, that is, a foam-filled shell. - In figure 
3 a general cross section for such a cylinder is illustrated. 
Because the thickness of the shell wall is small compared 
with the overall dimensions of the cross section, the stress 
in the wall can be assumed to be uniformly distributed over 
the thickness. This stress is equal to the normal derivative 
and is given by 

Ajl a I 

(ID 


btj>l 


1 &F 

t 


where represents the stress in the shell wall. With the 
use of equation (11), equations (4), (5), and (7) can now be 
written as 


and 




( 12 ) 


dr 


</n 

Gi & 

Cj ^1 t 

Cj t 




(13) 


Let the curve G in figure 3 be the middle surface of the 
shell wall. If the shell wall is assumed to be concentrated 
at its middle surface, then C can be thought of as represent- 
ing the interface, middle surface, and outer boundary of the 
wall. Consider R as being the region bounded by G. The 
problem can now be formulated as follows. Find a function 
<f> satisfying the equation . 


vV=— 2(? 2 0 

in R and the equation 

i>4> 
bv 



(14) 


(15) 



Figure 3. — Cross section of a composite thin-walled cylindrical shell. 


along C. The moment on the cross section is equal to 
twice the volume under the 4> diagram, or 

M=2jj 4>dxdy (16) 

ENERGY APPROACH 

Approximate solutions for the torsion of foam-filled shells 
can be obtained by the energy method. The complementary 
energy for such a body is the sum of the stress energy of the 
core, the stress energy of the shell wall, and the negative 
of the work done by the external moment acting through 
the angle of twist. For the composite cross section shown 
in figure 3 the complementary energy per unit length is 

U= + T v*)dx dy+^JVwan)^ ds—Me (17) 

where t„ and r„ represent the shear stresses in the core 
and r wm jj represents the shear stress in the shell wall. In 
terms of the stress function <f>, U becomes 

u =^JJ R (4>r 2 +^-iGM)dx ^+2^J/| c y (18) 

where the subscripts on $ denote the partial derivative with 
respect to the indicated variable. 

When the variation of Z7is equated to zero and integrations 
by parts are carried out, the following equation is obtained: 

O 3 SU=0^Jj^+^H ds-jj^+^G^ dx dy 

(19) 

It is seen that, if 5<f> is arbitrary in R and along C, equations 

(14) and (15) must be valid. 

The torsion of a foam-filled shell is analogous to the 
problem of the deflection of a membrane stretched over the 
region R, subjected to lateral pressure, and supported along 
the curve G by infinitely many springs which are constrained 
to distort only in the direction normal to the plane of R. 
Some discussion of approximate solutions of problems of this 
type is given in reference 4. 

SOLUTION OF SPECIFIC PROBLEMS 

In this part of the report an exact solution for a rectangular 
cross section is obtained by satisfying equations (14) and 

(15) . For the slender triangular and slender diamond cross 
sections, exact solutions do not appear feasible, and ap- 
proximate solutions are obtained instead by using equation 
(18). Two approximate procedures are utilized in each 
case; the first is the Rayleigh-Bitz method and the second is 
a more general variational procedure, herein called the 
“variational method.” This latter procedure is applied, 
for the most part, in cases where the Bayleigh-Ritz method 
becomes cumbersome. These two approximate methods 
may be applicable to other sections of. practical interest for 
which solutions are not available. 
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RECTANGULAR CROSS SECTION 

Tlie notation for a rectangular cross section is shown in 
figure 4 (a). The thickness of the shell wall is assumed to 
be constant. The problem is to find a function <f> which 
satisfies equation (14) within the rectangle and the following 
boundary conditions: 

iHf? (m 

The function <f> is, of course, symmetric about both the X- 
and T-axes. The symmetry condition about the X-axis is 
satisfied when <j> is taken in the form 

COS X 4/ ' (21) 

»=0 


Y 



(a) Rectangular cross section. 

(b) Triangular cross seotion. 

(c) Diamond cross section. 
Figure 4. — Notation used in analysis. 


where the functions cos X^y form an orthogonal set. The 
right-hand side of equation (14) can be expanded in a Fourier 
series of the functions cos X„y in the interval —b^y&b, and 
this expansion yields 

.. 40 s 0 sin X*£> • . 

' V s 0=— 2-1 7 , ■ \ I . i cos X*y (22) 

o X„o-f- sm \J> cos X*o 

Substituting the assumed solution (eq. (21)) into equation 
(22) and equating coefficients of like terms gives the following 
ordinary differential equation. for h x (x): 


h n (x)"-\SK(x)=- 


4 Ojd sin \ H b 
X,6-f-ain X„& cos X,6 


The solution to equation (23) is 
h n (x) =k ln sinh \ n x+k iH cosh X„a;+ 


4020 sin X a 6 

^» a (X»6+sin \ n b cos X„6) 

(24) 


The constants k ln vanish because of symmetry. The con- 
stants k lx and the eigenvalues X* can be found from the 
boundary condition^ at x=a and y=b, respectively. Con- 
sider first the condition at y=b. The substitution of equation 
(21) into equation (20b) yields 

X B sin X„&=i^ cos \ K b (26) 

Therefore, the eigenvalues are given by 

tan X b 6=^ (26) 

where 


At x=a, the substitution of equation (21) into equation (20a) 
yields 

7 . • -» . 0* i . . 40j0 sin \ n b "1 

fe.X. smh X.«— Jy |fe, cosh X.«+ V(Xji+BillX<icosX<t) j 

(27) 

Therefore, 

, —^KGjdb* sin X,6 

Ss— X„*& S (X*& sinh X^a+X cosh X*a)(X B 6-l-sin \ n b cos X„6) 

(28) 

Consequently, the stress function is 

^ 40j/db 3 sin X,5 cos X»y 

*-£* \, i b i (k H b + sin X*6 cos X*6) X 


K cosh X»a: 

X,Z> sinh \a-\-K cosh X 


s) (29) 


The moment on the cross section is given by the formula 
M= 8 f f <j>dxdy (30) 
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The torsional stiffness can be expressed in terms of either O x 
or 0 2 

-t=G\J\=QiJ* (31) 


The torsional stiffness constants J\ and J, are obtained upon 
substituting equation (29) into equation (30) and carrying 
out the indicated integration. When the results are expressed 
in the form of equation (31), it is seen that 


where 


J\=Act^tKk' 

J,=2ct 0 3 A 


A =y^ . &m 2 \ n b 

jpo X* 3 X 3 (X n X+sin \ K b cos \»b) 


[l-- 

a 


K sinh X,g "1 

X«6(X,6 sinh \ n a+K cosh X^cOJ 


The series A converges very rapidly. For various values of 
K, the eigenvalues are easily located from the intersections 
of the hyperbola K/\ n b and the curves tan X,6 as indicated 
in figure 5. Plots of <7, and J, against the cross-section aspect 
ratio U/c for various values of K are presented in figure 6. 

As the stiffness of the shell-wall material vanishes, K 

approaches infinity and X„6 approaches It can be 

shown that the limiting value of J, is the solution obtained 
by the theory of elasticity for the torsion of a homogeneous 
rectangular cross section as given in reference 3 (p. 278). 

As the stiffness of the core material vanishes, K approaches 
zero and the solution should approach that of the Bredt 
theory for the torsion of hollow thin-walled shells with free 
warping. As K approaches zero, \J) approaches mr. It is 
obvious, then, that all the terms in A vanish for which w]>0. 
Investigation of the term for which n = 0 shows that J x does 


approach the result given by the Bredt theory which is 

v T 4A 2 2 ct 0 H 

lim ° 

-sr-o 7, , 


fj K) 


where A 0 is the area enclosed by the median line of the shell 
wall. 

SLENDER TRIANGULAR CROSS SECTION 

If the energy approach is used, two approximate solutions 
are obtained for a cross section in the shape of a slender 
isosceles triangle with a constant-thickness shell wall as 
illustrated in figure 4 (b). One solution is obtained by the 
Rayleigh-Ritz method and another solution is found by 
utilizing the calculus of variations and the boundary-layer 
technique of reference 5. 

Eayleigh-Kitz method. — In terms of nondimensional co- 
ordinates, the complementary energy (eq. (18)) can be 
written for the triangle as follows: 

17=0,11 

= f f W+V-4 f VI d£+ 

Jo Jo £ £jo|ij=f 

-ST f V dti (35) 
Jo |£~1 


where 




Note that the equation of the sloping side of the triangle in 
nondimensional coordinates is simply 17= £. 

The stress function <f> must be an even function through 
the thickness, and for slender sections it is usually sufficient 
to assume a parabolic variation in the thickness direction. 
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For thick sections, however, it may be necessary to include 
additional terms in the thickness direction to get satisfactory 
accuracy. Suppose now that 4> is assumed to be a poly- 
nomial 


and the following boundary conditions are obtained: 

($/'+ (42a) 


n^Q 


(36) 




(42b) 


When equation (36) is substituted into equation (35) and U 
is minimized with respect to the parameters am and j9, the 
following r+2 equations result: 


(ef+Kf)^=0 (42c) 

(<g'+'.Kg)t-i=0 (42d) 


Z>U 

fern 


=0 


n«0 \ 


2 mn ^ 


Kp. 


m+n m+n - hi 


2 


to +2 


GjV 


W 

2>J8 


=0 


r (37) 


=S“-(55+T)+<f + f+lH W ’ 

where 

Solution of the system of simultaneous equations (37) 
yields a, and /3. From equation (16), the moment is 


M =<4SiT2+i) 


(38) 


The stiffness can be written in the form of equation (31) and 
Ji and J t are easily calculated. The results are 



1 «» 
7J.+2 Q$b 2 



I _ c< o 3 



1 ■ P_ 

n+2 G s 6b 3 ' r 2G 3 eb 3 


-) 


f (39) 


When K is large compared with unity a large number of 
equations may be required for reasonable accuracy. An 
approach which avoids this difficulty is developed in the 
succeeding section. 

Variational method. — Instead of assuming for 4> the poly- 
nomial of equation (36), suppose arbitrary functions of £ are 
allowed to remain and <j> is taken to be of the form 

4> ==/(£) +/i;(£) (40) 

When equation (40) is substituted into equation (35) and 
the variation of TJ with respect to admissible variations in 
/ and g is equated to zero, two simultaneous ordinary differ- 
ential equations for the functions / and g are obtained as 
follows: 


y+f (ftr-w+ft)= -2 ^ 


k (i s /')'+f (.m'-Kum+te)— ojb 2 ? 


( 41 ) 


where the primes denote derivatives with respect to £. 

The differential equations (41) are linear with variable 
coefficients, and it appears to be a difficult task to find an 
exact solution to the system. For the case of slender cross 
sections, however, an approximate solution is possible by 
utilizing the “boundary layer” technique discussed in refer- 
ence 5. Notice that the differentiated terms in equations 
(41) are multiplied by £, a quantity which for slender cross 
sections is small compared with unity. Differential equa- 
tions having the most highly differentiated terms multiplied 
by a small quantity are characteristic of the type of bound- 
ary-layer problems considered in reference 5. 

Suppose, initially, that/ and g are slowly varying functions 
throughout the region 0 ^ £ A 1 . The term “slowly varying” 
is intended to mean that the maximum values of the func- 
tions/ and g and their derivatives which appear in equations 
(41) are of the same order of magnitude. Then, as long as 
K is at least of the order of unity, the terms in equations 
(41) which contain e 2 have little influence on the solution. 
Consequently, a good approximation to a particular solution 
is obtained by ignoring the terms in equations (41) which 
contain £. 

When this procedure is carried out, it is found that the 
approximate particular solution satisfies the boundary con- 
ditions at £=0 but does not satisfy the boundary conditions 
at £=1. It can be concluded that the required solution is 
such that the functions / and g are not slowly varying every- 
where in the region 0^£^1. Somewhere the derivatives 
which appear in equations (41) must take on values which 
are of the order of e~ 2 so that the terms containing £ can 
have an appreciable influence on the solution. 

It is assumed that the region where the derivatives of 
/ and g have values of the order of e~ 2 is confined to a so-called 
boundary layer in the neighborhood of £=1. On the basis 
of this assumption, the particular solution alone is a good 
approximation to the exact solution away from £=1. Then, 
by focusing attention on the boundary layer close to £=1, 
it is possible to obtain an approximate homogeneous solution 
to equations (41) which modifies the particular solution in 
such a manner that the boundary conditions at £=1 can be 
satisfied. 

It is convenient to get a particular solution as a power 
series in e instead of ignoring completely the £ terms in 
equations (41). Assume that a solution can be expressed 
in the form 


/p(5)=/po(f)+«fpi (£)+«/«(£)+ • • • 
5e(£) —giti (£) + tgpi (£) + (£) + • • • 


(43) 
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where the subscript P denotes a particular solution. When 
equations (43) are substituted into equations (41) and 
coefficients of like powers of e are equated, pairs of simul- 
taneous equations are obtained for the coefficients in the 
power series. For example, when the coefficients of the 
zeroeth power of e are equated there results 

^m(/po+£ 2 £po)=2(? 2 0& 2 £ 



When coefficients of the first power of e are equated the 
result is 


fn'ptfgpi — 0 

•^M(/pi+£ 2 £n)4~g fepi— 0 


(44b) 


When coefficients of the second power of e are equated the 
following equations result: 


(£/po')' +Tj (£ 3 gi’o / ) / — -K)*(/p*+f s £p2) — 0 
2 (S 3 : /poO'+g (l?gp o')' — — g f 3 £p2=0 


H44c) 


Similar equations are obtained when coefficients of higher 
powers of e are equated. 

Solution of equations (44) results in the following expres- 
sions for j P and g P : 


f P -6jV [_J-+¥+ S (gyH A+ 1 **) + " '] „ 

£p =-Gr^ (1+£i+ . . .) 


(46) 


Coefficients of the odd powers of e vanish. This particular 
solution satisfies the boundary conditions at £=0 but not 
at £=1. 

A homogeneous solution can be obtained which modifies 
the particular solution in the vicinity of {=1 in such a way 
that the boundary conditions at £=1 can be satisfied. In 
order to determine the homogeneous solution it is convenient 
to put equations (41) into a form in which, in the neighbor- 
hood of £=1, the terms containing derivatives are of the 
same order of magnitude as the remaining terms. Such a 
conversion is provided by the coordinate transformation 

£=1+4 (46) 


the following equations are obtained: 

F(i+4)/ , ] , 44[(i+4)y] / -^[/+(i+4) 2 ^]=o' 


f[(i+©yt+f[(i+©Y]'-i^[(i+©’/+ }■ W 

(i+4) 4 d-| (i+4) 3 s=o 


where the primes now denote differentiation with respect 
to £. 

The homogeneous solution can bo expanded in powers of e 
fs(£)=fm(£) + «/?i(D +e%2© + • • •") 

r (48) 

&r©=&o© + e&i© + e 2 fo2© + * • *J 


where the subscript H denotes a homogeneous solution. If 
these expressions are substituted into equations (47) and the 
coefficients of like powers of e are equated, systems of 
ordinary differential equations are obtained. For the terms 
/so and gm, a set of homogeneous equations is obtained: 


ja o // + ^gao " — -E"m (/ao+gflo) — 0 
g/m> // + (fan~\~gm)~ gg//o~0 


(40) 


For the terms f m and g, n the following nonhomogeneous 
equations are obtained: 


/t 2 // + qSbi" — Kp- (fm +&n) — — (1/voO ' — Qga o') ' + 

2Kplgm 

^Sai" — Kp(ftnPgm) — g gai — — (ijfcoO ' — 
Qgm') ' + (/so + 2 gaa) + ^kgno 


(60) 


Additional sets of nonhomogeneous equations would result 
for the coefficients of higher order terms. 

It is found that neglecting terms of the order of e in equa- 
tions (48) is equivalent to neglecting terms only of the order 
of e* in the final result for the torsional stiffness constants. 
Therefore, a final result which includes all terms linear in 
e can be obtained by solving only equations (49) and dropping 
all higher order terms in the homogeneous solution. Solu- 
tions to equations (49) are of the form 


When the transformation (eq. (46)) is introduced into 
equations (41) and the right-hand sides are set equal to zero, 


f m =Ae* 
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Substitution of these functions into equations (49) yields 


AiV-Krf+B^-K^O 

a($-k,)+b(?-k^)-o 


y (si) 


.and where terms of the order of e 2 have been dropped. It 
is seen now that in calculating the arbitrary constants A n , 
terms of the order of f may be neglected. When equations 
(54) are substituted into the boundary conditions (eqs. 
(42c) and (42d)) the following expressions for A n are obtained 
after dropping terms of the order of e: 


This system of equations has a solution only if the deter- 
minant of the coefficients vanishes. When the determinant 
is set equal to zero a biquadratic equation for X is obtained. 
The solution to the biquadratic equation is 


^ ^ 15-\-QKn i -V(15+6.K/x) 2 — 60Kp 


(52) 


Q+E^-2K 

(Xi-b-Zr)V(15-[-6l£)i) 2 — 60 Kn 

(~A-Ky^-2K 

(X2 +^)V( 15 + 6.K71) 2 — 60 -K/i 


(&t 


The homogeneous solution must vanish as the absolute 
magnitude of £ increases, and since f is always negative or 
zero, only the two positive roots for X are required. 

From the first of equations (51) the B coefficients can be 
written in terms of the A coefficients 


B =- 3 


x *~ K m a 
X 2 -3 Kp 


(53) 


When the £ coordinate is transformed to the £ coordinate and 
terms of the order of e 2 are dropped, an approximate general 
solution to equations (41) is obtained 


K V 


f—^2A x e 'e* 

S=-2 3 A x e~e^-Gjb'- 

n~l X n — Oit/i 


(54) 


This solution satisfies the boundary conditions at £=0. 
The boundary conditions at £=1 are used to determine 
Ai and A 2 . 

When equations (54) are used, the stress function <p 
becomes 


In the limiting case where K approaches zero the boundary- 
layer technique becomes invalid. Therefore, the solution 
cannot be expected to approach the proper result for a hollow 
thin-walled shell. When K approaches infini ty, an approx- 
imate solution for a solid cross section is obtained; and T 
is given by 


lim T=l— 4 

K—*CO 



€ 


(59) 


SLENDER DIAMOND CROSS SECTION 

Bayleigh-Bitz method. — For the slender diamond cross 
section the notation is shown in figure 4 (c). The comple- 
mentary energy becomes 


dv d *+ 2 tICH (eo) 

Substitution of the polynomial (eq. (36)) into equation (60) 
and minimization with respect to a m and jS yields the r+2 
equations 





1 — 


3 


XjLZjfr 

Xn 2 — 3 Kfl 





(55) 


where 



and where terms of the order of e 2 are dropped. The moment 
is calculated by substituting 4> into equation (16) and per- 
forming the integration. The results for the torsional 
stiffness constants J, and J 2 are 


where 


Jr- 


ct 0 3 t 
= 12 

j-cjl r 
12 r 


2KT 


r=l+ry — b6 An 

•R-ft n«=I 



(56) 


(57) 


i)U . f , mn 

o&ot n-o V m+n 


m-fw+l 



#M 

m+3 


§=0 =2>„ 

Op n-0 



2 

m+2 


GsdbA 


y 


J 


(61) 

where O^m^r. 

The moment and the torsional stiffness? constants J x and 
J 2 are found by utilizing equation (16). The expressions 
for the torsional stiffness constants turn out to be precisely 
the same as for the triangular cross section given in equations 
(39). Of course, for the diamond cross section, a„ and are 
obtained from equations (61). 

Variational method. — Through the use of the calculus of 
variations and expression (40) for <f>, the differential equations 
(41) are found to be valid also for the diamond cross section. 
The boundary equations at £=0 (eqs. (42a) and (42b)) also 
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hold for the diamond cross section. The boundary condi- 
tions at £=1, however, are now given by 


In the case where K approaches infinity an approximate 
solution for a solid section is obtained, and T* is given by 


/'U=ol 
£'li-i=o J 


where the primes denote differentiation with respect to £. 

The boundary-layer technique yields a general solution of 
the same form as that obtained for the triangle. However, 
now it is found that neglecting terms of the order of e in the 
homogeneous solution (eqs. (48)) is equivalent to neglecting 
terms only of the order of £ in the final expressions for the 
torsional stiffness constants. Consequently, for the diamond 
cross section a result which includes all terms of the order 
of e 2 can be obtained by solving equations (49) and dropping 
all higher order terms in the homogeneous solution. It is 
consistent, now, to keep the e term in the particular solution. 

The approximate general solution to equations (41) for 
the diamond cross section becomes 






y (63) 


The arbitrary constants A x and A 2 are determined by sub- 
stituting equations (63) into equations (62). 

The stress function becomes 


If this expression for <f> is substituted into equation (16) and 
integrated, the moment can be calculated. The torsional 
stiffness constants become 


where 


j- c -h^2Kr* 
12 1 


(65) 


r*=i 




( 66 ) 


When equations (63) are substituted into the boundary 
conditions (eqs. (62)) and terms of the order of e 2 are 
neglected, the following expressions are obtained for the 
arbitrary constants: 

-q- 1+Kn 2X 2 2 -f 3jST/i 

Xj -\/(l 5 + 6 Kfif — 60K/J. 

" (67) 

-j__ 1 +Kn 2X,»+3 Ky 

■K/* X j V(15+ 6 Kitf — fjOKp __ 


lim r*=l-3e 2 (68) 

j5T— 


RESULTS AND DISCUSSION 

PRESENTATION OF RESULTS 

The results of the calculations outlined in this report are 
presented in figures 6, 7, and 8. These figures show plots 
of the torsional stiffness constants against cross-section 
aspect ratio for various values of the parameter K. Tho 
torsional st iffn ess can be expressed either in terms of the 
shear modulus of the core material or the shear modulus 
of the shell-wall material. The torsional stiffness constant- 
associated with the shear modulus of the core material Jj 
is plotted in parts (a) of figures 6, 7, and 8, whereas the 
constant associated with the shear modulus of the shell-wall 
material J x is plotted in. parts (b) of these figures. 

In figure 6 are shown the results of the exact solution of 
the differential equation (14) with the boundary conditions 
(20) for a rectangular cross section. In figures 7 and 8 are 
shown the results of the approximate solutions (the Rayleigh- 
Ritz method and the variational method in conjunction 
with the boundary-layer technique) for the triangular and 
diamond cross sections. A five ' parameter polynomial 
was used in the Rayleigh-Ritz method. 

accuracy of approximate methods 

Solutions by the Rayleigh-Ritz and variational methods 
also were obtained for the rectangular cross section. A 
comparison of these results with the exact solution provides 
an indication of the accuracy of the approximate methods. 

A polynomial with only three parameters was used for 
the Rayleigh-Ritz method in this comparison. The results 
showed that for all aspect ratios and for K less than about 
unity the stiffness given by Rayleigh-Ritz method is loss 
than 3 percent lower than the exact stiffness. It is believed 
that the five-parameter Rayleigh-Ritz method used for tho 
triangular and diamond cross sections should yield slightly 
more accurate results. Of course the accuracy of tho 
Rayleigh-Ritz method can be improved for large values of 
K. by including more terms in the polynomial for <j>. How- 
ever, the number of simultaneous equations which must bo 
solved increases with the number of unknown parameters. 

The stiffness calculated by the variational method for 
the rectangular cross section was less than 1 percent in error 
for values of K greater than about unity and tho aspect 
ratio t 0 fc less than about K- The boundary-layer technique 
yielded a slightly more accurate solution to the differential 
equations obtained for the rectangular cross section than for 
equations (41) which arise for the triangular’ and diamond 
cross sections. Thus the results of the variational method 
for the triangular and diamond cross sections are probably 
not quite as accurate as for the rectangular cross section. 
It appears that for slender cross sections (small values of 
fo/c) the variational method is more accurate for largo 
values of K and the Rayleigh-Ritz method is more accurate 
for small values of K. . 




Royleigh-Ritz method 

Vorlotionol method 















TOBSION OF THIN-WALLED SHELLS WITH REINFORCING CORES 


In figures 7 and 8 the solid curves axe results of the five- 
parameter Rayleigh-Ritz method, and the dashed curves 
show results of the variational method. The Rayleigh- 
Ritz method certainly leads to a lower bound for the torsional 
stiffness. The variational approach also leads to a lower 
bound provided the boundary-layer technique yields a 
sufficiently accurate solution to equations (41). From the 
results on the rectangular cross section it appears reasonable 
to assume that for IsTe^I the variational method probably 
gives a lower bound. Therefore, for any particular case 
where 1, the higher of the two values of torsional stiff- 
ness calculated by the two approximate methods is the more 
accurate. 

BEMAKKS ON AN ELEMENTARY CALCULATION 

It is conceivable that a first approximation to the stiff- 
ness of a foam-filled shell could be made simply by adding 
the torsional stiffness of the core alone to that of the shell 
alone and neglecting the stiffening effect which results from 
bonding the core and shell together. Calculations were 
made by this elementary procedure and the results are shown 

in figure 9. The ratio is the ratio of J l as calculated 

by the elementary procedure (neglecting the bonding effect) 
to Ji as calculated by the theory of this report. It is seen 
that for the rectangular cross section the result of the ele- 
mentary calculation is never more than 5 percent lower 
than the exact solution. For the triangular and diamond 
cross sections, however, the elementary procedure does 
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not lead to such good results, and the discrepancy can be 
as much as 25 percent. 

CONCLUDING REMARKS 

An approximate boundary-value problem is set up for the 
Saint-Venant torsion of cylindrical thin-walled shells bonded 
to a core which fills the interior of the shell and which is 
made of a material different from that of the shell wall. 
Solutions for the torsional stiffness are obtained for three 
particular cross-sectional shapes — an exact solution to the 
boundary-value problem for rectangular cross sections and 
approximate solutions for slender triangular and slender 
diamond cross sections. The approximate solutions are ob- 
tained by the use of two energy procedures. These methods 
may be applicable to other cross sections of practical interest. 
The choice of approximate method for any particular prob- 
lem depends on the range of parameters involved. 

The stiffness obtained by the simple procedure of adding 
together the individual stiffnesses of the core and the hollow 
shell (neglecting the effect of the bond) yields results less 
than five percent low for rectangular cross sections. For 
slender triangular and diamond cross sections this elementary 
approximation is generally not so good and in certain cases 
it yields results which are considerably low. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 7, 1956. 
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(c) Diamond cross section. 

Fiqubb 9.— Comparison between the results for the torsional stiffness constant J i obtained from an elementary calculation and the present theory. 




